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Abstract 

Particle identification (PID) capabilities are studied by using the Time Projection 
Chamber (TPC) and a Time-Of-Flight (TOF) detector together at STAR. The 
identification capability of charged hadrons is greatly extended compared with that 
achieved by TPC and TOF separately. Particle spectra from p+p, d+Au collisions 
at y/s NN =200 GeV and Au+Au collisions at ^s NN =62.4 GeV are used to develop 
the methods. The transverse momentum (pr) ranges of tt, and p(p) identification 
are from ~ 0.3 GeV/c to ~ 10 GeV/c. The high pt reach is limited by statistics 
in current data sets. An important conceptual advance was developed to identify 
electrons by using a combination of dE/dx in TPC and velocity information from 
the TOF detectors, which is important for future low-mass dilepton program at 
STAR. 
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1 Introduction 



One of the goals of the relativistic heavy ion program at RHIC is to study 
quantum chromodynamics (QCD) at extreme condition [1]. A unique strength 
of the solenoidal tracker at RHIC (STAR) [2] is its large, uniform acceptance 
capable of measuring and identifying a substantial fraction of the particles 
produced in heavy ion collisions. Detectors relevant to the study presented 
in this article are the Time Projection Chamber (TPC) [3], and a proposed 
barrel time-of-flight (TOF) [5]. For stable charged hadrons, the TPC provides 
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tt/K (7r+K/p) identification to px — 0.7 (1.1) GeV/c by the ionization energy- 
loss (dE/dx) as usually been quoted and presented in the physics analyses [3]. 
However, direct particle identification (PID) capability for stable hadrons at 
intermediate/high p T is important for the study of collective flow and strong 
early-stage interaction in the dense medium formed in relativistic heavy ion 
collisions [1]. STAR PID capability can be further enhanced by the proposed 
TOF. A TOF system with a time resolution of < 100 ps at STAR is able to 
identify n/K (n+K/p) to pt — 1.6 (3.0) GeV/c, as demonstrated in Fig. 1 (see 
also [6] [7]). In addition, with relativistic rise of dE/dx from charged hadrons 
traversing the TPC at intermediate/high p T (>3 GeV/c) and diminished yields 
of electrons and kaons at this px range, we can identify pions and protons up 
to very high px (— 10 GeV/c) in p+p, p+A and A+A collisions at RHIC. An 
important conceptual advance was developed to identify electrons by using a 
combination of dE/dx in the TPC and velocity information from the TOF 
detectors. This has been used to measure charm yield via its semileptonic 
decay [8]. Electron identification and hadron rejection power will be discussed 
in detail. This provides a basic tool for the future dilepton measurements with 
the azimuthal 2n coverage of a barrel time-of-flight system. The proposed 
dilepton measurements will provide a penetrating probe into the new state of 
dense matter produced in central heavy ion collisions at RHIC, since leptons 
do not participate in strong interactions occurring during hadronization and 
freeze-out. 



2 Experiment Setup 



At STAR[2], the main tracking device is a TPC, covering full azimuthal an- 
gle and ±1.5 units of pseudo-rapidity. A dE/dx resolution of ~ 8% can be 
achieved by requiring the tracks of charged particles to have at least 20 out 
of a maximum of 45 hits in the TPC. Detailed descriptions of the TPC and 
its electronics system have been presented in [3] [4]. One tray of prototype 
time-of-flight detector based on multi-gap resistive plate chambers [9] (TOFr) 
was installed in STAR in 2003. It covers 1/60 in azimuth and — 1 < rj < in 
pseudo-rapidity at the outer radius of the TPC, 220 cm from the interaction 
point. Two identical pseudo-vertex position detectors (pVPD) were installed 
to record the start time for the TOFr, each 5.4 m away from the TPC center 
along the beam line. Each pVPD covers ~ 19% of the total solid angle in 
4.4 < \r)\ < 4.9. More information of the TOFr and the pVPD can be found 
in [10] [11]. The data used in this study were collected from Au+Au collisions 
at y^SNN — 62.4 GeV in 2004, and p+p and d+Au collisions at ^^snn — 200 
GeV in 2003 at RHIC. The time resolution of the TOF system is ~ 105ps 
(Au+Au) and ~ 120ps (d+Au) respectively, which include pVPD's contribu- 
tion of ~ 55ps (Au+Au) and ~ 85ps (d+Au). 
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3 Stable Hadron Identification 

3. 1 PID at intermediate/high p T by TPC 



TOFr PIP (62 GeV Au+Au run) 




3.5 

p (GeV/c) 



Fig. 1. 1/(3 vs. momentum for pions, kaons and (anti-)protons from TOFr at 62.4 
GeV Au+Au collisions. The separation between pions and kaons ( (anti-)protons) 
is achieved to px ~ 1-6 (3.0) GeV/c. 
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Fig. 2. Distribution of logio(dE jdx) as a function of logio(p) for electrons, pions, 
Kaons and (anti-)protons. The units of dE/dx and momentum (p) are keV/cm and 
GeV/c, respectively. The color bands denote within ±la the dE/dx resolution. 170 
means Bichsel's prediction for 30% truncated dE/dx mean. 

At 3 < Pt^IO GeV/c, there is a difference of about 15% in the dE/dx between 
pions and kaons due to the pion relativistic rise of the ionization energy loss. 
The difference between that of pions and (anti-)protons is even larger. This 
allows us to identify pions from other hadrons at this px range by the TPC 
alone at 2a level. The dE/dx resolution is ~ 8%, as demonstrated in Fig. 2. 
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Fig. 3. dE/dx distribution normalized by pion dE/dx at 4 < px < 4.5 GeV/c and 
\r/\ < 0.5, and shifted by ±5 for positive and negative charged particles, respectively. 
The distribution is from minimum-bias Au+Au collisions at ^/s NN =62.4 GeV. 
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Fig. 4. The relative dE/dx peak position of K-tt (squares) and p-7r (cir- 
cles) in unit of standard resolution width (a n ) of pion dE/dx as function 
of pt- The triangles are the peak positions of the dE/dx distribution of 
h + - h~ = (p - p) + (K + - K~) + (7T+ - 7r~) ~ ip-p). 

Shown in Fig. 3 is the no^ distribution for charged hadrons at 4 < px < 4.5 
GeV/c and \r]\ < 0.5, where na n is the normalized dE/dx of pions. The nor- 
malized dE/dx is denned by na x = log((dE / dx)y / B x ) / ox, m which X, Y 
can be e ± ,7T ± ,K ± or p{p). B x is the expected mean dE/dx of a particle X, 
and ox is the dE/dx resolution of TPC The na^ distribution is a normal 
Gaussian distribution with an ideal calibration. The na n of positive and neg- 
ative charged hadrons are displaced by +5 and -5 respectively in Fig. 3. Fig. 4 
shows the px dependence of na^ and na^ relative to ncxj, as predicted by the 
Bichsel function for the energy loss in thin layers of P10 [3]. From Fig. 3 and 



4 



measurement from TOF [7], the yield difference between positive and negative 
inclusive charged hadrons is approximately that of proton and anti-proton ( 
h + — h~ = (p — p) + (K + — K~) + — 7r~) ~ (p — p), kaons' contribution 
is small). Therefore, the peak positions of dE/dx distribution from h + — h~ 
should represent well that of protons (since na% + ~ h ~ na^^ = na^). 
Indeed, the calibrated Bichsel function for protons matches the dE/dx peak 
position of (h + — h~), as shown in Fig. 4. In addition, as shown later, the 
dE/dx difference of protons and pions in the momentum range where PID 
selection is possible by TOF is also found to be consistent with the Bichsel 
function. These crosschecks confirm that the Bichsel function can be used to 
constrain the relative dE/dx position between kaons, (anti-)protons and pions. 
To extract pion yield, we performed a six Gaussian fit to the dE/dx distri- 
butions of positive and negative hadrons simultaneously as shown in Fig. 3. 
The na^ — na^ and na^ — na^ are fixed in a six- Gaussian fit, where the six 
Gaussians are for 7r ± , K ± and p{p) at a given p T bin. The sigma of the six 
Gaussians are chosen to be the same. The pion yields extracted from the fit 
can be found in [12]. As shown by Fig. 2 and discussed in [12], the PID of 
pions at ~ 3 < px <~7-8 GeV/c is obtained with this method. The pr reach 
is limited by statistics and not by PID capability in this study. 




Fig. 5. h /h + vs dE/dx distribution normalized by pion dE/dx at 3.5 < pt < 4.0 
GeV/c and \r)\ < 0.5. The distribution is from minimum-bias Au+Au collisions at 
= 62.4 GeV. 

Fig. 4 shows that the dE/dx separation between kaon and proton is less than 
one <7 at pr between 3 and 5 GeV/c and larger at pr around 10 GeV/c. How- 
ever, there are a few methods we can use to identify protons and cross check 
the contaminations from kaons. Fig. 5 shows the ratio of negative hadrons 
over positive hadrons as a function of dE/dx in unit of na n at 3.5 < pt < 4.0 
GeV/c. Since the energy loss of particles in TPC is independent of its charge 
sign, the dependence of h~ /h + on na w is due to different particle composition 
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Fig. 7. m 2 distribution without dE/dx 



Fig. 6. Particle distribution as 

a function of nov and m 2 for , 

. _ xr/ . . cut ( black solid line), with nov > U 

2 < < 4 GeV/c. The pion, kaon ' , „ 



and (anti-)proton peaks can be 
clearly seen. 



cut (red dashed line), and na p < 
cut (blue dotted line). 



and due to the dE/dx separation between pion, kaon and proton. Fig. 5 shows 
two plateaus. One from p/p and the other from tt~ /t[ + . We will be able to se- 
lect protons with reasonable purity by requiring na^ < —2.5 in this particular 
Pt bin. We can then use an independent measurement of Kg from VO method 
to constrain the charged kaon yields in the six-Gaussian fit or to estimate the 
charged kaon contaminations described above to obtain reliable proton yield. 
Detailed studies and measurements are underway. 



3.2 Hadron PID at Intermediate p T by TPC+TOF 



Fig. 1 and Fig. 2 show that TOF alone is not able to separate between pions 
and kaons and TPC alone is not able to separate pions, and protons at the 
Pt range roughly between 2.0 and 4.0 GeV/c. However, the dependence of 
1/(3 on pt from TOF and that of dE/dx on pt are different in this pt range 
and therefore, by combining these two, we will be able to extend our PID 
capability. Fig. 6 shows the hadron distribution as a function of na n and mass 
square (m 2 ), with m 2 = p 2 ((tTOF * c/l) 2 — 1), where p is the momentum, t T oF 
is the time of flight, c is the speed of light in vacuum, and I is the flight path 
length of the particle. The pion, kaon and (anti-)proton peak can be clearly 
seen in Fig. 6. The solid line in Fig. 7 is the projection of Fig. 6 to the m 2 
direction at 3 < pt < 4 GeV/c. At this pt range, the pion and kaon bands 
are merged together, and cannot be clearly separated from the (anti-)proton 
band with the TOF. However, if we require na n > and then plot the m 2 
distribution again (dashed line in Fig. 7), the kaon and (anti-)proton bands are 
greatly suppressed and a clear pion signal is observed. This is due to the rather 
large difference between na^ and na^ or n<j p p ^ in this pt range, as shown in 
Fig. 4. Similarly, the pion and kaon bands are suppressed significantly with 
respect to the (anti-)proton band when na p (p) < is applied. This helps us 
get a cleaner (anti-)proton signal. Therefore, the combination of the TPC and 
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TOF enhances the particle identification at 2 < p? < 4 GeV/c where neither 
TPC nor TOF can well separate the hadrons. 

The invariant yields of pions, kaons and (anti-)protons were calculated for 62.4 
GeV Au+Au collisions. The results were shown elsewhere in [7]. The pions and 
(anti-)protons were identified at 0.2 < pr <~ 5 GeV/c, and the kaons were 
identified at 0.2 < px <~ 3 GeV/c. The PID px reach is due to low statistics 
with the small acceptance of TOFr. 



3.3 Cross-check at low /intermediate pt 



pT range: (1.400000,1.600000) 



Entries 26112 
Mean 0.3185 
RMS 0.3995 
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Fig. 8. to 2 distribution from TOFr for 
1.4 < pt < 1-6 GeV/c. Arrow shows 
the cut at the pion mass in order to 
obtain a clean pion sample for dE / dx 
studies. 




Fig. 9. na n distribution using the 
pion sample resulting from the cut de- 
scribed in previous figure. The solid 
lines are fits to a Gaussian function. 



Since dE/dx plays a key role in our methods discussed above, dE/dx calibra- 
tion has to be carefully studied. Again, with PID capability at overlapping px 
range, we can use identified particles (pions, protons etc.) from TOFr to cross- 
check the characteristics of na{n, K, pip)) distribution at low/intermediate px- 
Fig. 8 demonstrates an example of choosing a pion sample by cutting on m 2 . 
We can then plot the na^ distribution of this sample and fit it with a Gaus- 
sian function as in Fig. 9. The mean and sigma of the Gaussian are listed in 
Table 1 as a function of px- The parameters (mean and sigma) obtained from 
proton samples are also listed in Table 1. We note that the width of ntrj is 
systematically larger than that of pions, partly due to the different resolu- 
tion from different track length resulting in dispersion of na?. The values of 
these parameters for pions and (anti-)protons are almost independent of px 
and consistent with that obtained by pure dE/dx method (see Fig. 3). The 
last two columns in Table 1 show the measured < na^-ncr* > compared 
to Bichsel function calculation. The experimental results agree well with the 
expectations within errors. 
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Table 1 



na n Gaussian fit parameters of the pion and (anti-)proton sample obtained using 
m 2 cut. The differences between n<r%^ and na% are also listed and compared with 
Bichsel calculations. At pr > 3 GeV/c, < na^-na 7 ^ > are calculated via h + — h~ . 



Pt 
(GeV/c) 




p(p) 






mean 


sigma 


mean 


sigma 


measurement 


expectation 


0.9-1.0 


-0.185±0.006 


0.980±0.004 


4.727±0.014 


1.110±0.010 


4.912±0.020 


4.936 


1.0-1.2 


-0.171±0.005 


0.977±0.004 


3.176±0.021 


1.171±0.015 


3.347±0.026 


3.421 


1.2-1.4 


-0.150±0.007 


0.958±0.005 


1.562±0.013 


1.112±0.011 


1.712±0.020 


1.872 


1.4-1.6 


-0.180±0.010 


0.946±0.007 


0.464±0.015 


1.073±0.012 


0.645±0.025 


0.768 


1.6-1.8 


-0.195±0.013 


0.938±0.010 


-0.322±0.019 


1.021±0.014 


-0.127±0.032 


-0.036 


1.8-2.0 


-0.228±0.021 


0.923±0.015 


-0.834±0.023 


0.989±0.016 


-0.606±0.044 


-0.616 


2.0-2.5 


-0.171±0.029 


0.972±0.023 


-1.458±0.025 


1.034±0.019 


-1.287±0.054 


-1.174 


2.5-3.0 






-1.912±0.056 


1.161±0.048 


-1.91±0.01 


-1.72 


3.0-3.5 






-2.387±0.113 


1.111±0.121 


-2.40±0.01 


-2.06 


3.5-4.0 










-2.50±0.02 


-2.27 


4.0-4.5 










-2.52±0.04 


-2.35 


4.5-5.0 










-2.47±0.08 


-2.44 


5.0-6.0 










-2.47±0.15 


-2.47 


6.0-7.0 










-2.60±0.28 


-2.47 



3.4 Contamination in pion identification 



rdE/dx ir Identification, 4.0<p T <4.5 GeV/c 




Efficiency 

Fig. 10. Contamination from Kaons+(anti-)protons to pions with a dE/dx selection 
window vs pion efficiency from the same selection at 4 < pt < 4.5 GeV/c. Dashed 
line is for ir + and solid line is from 7r~. 

Good particle-by-particle identification is crucial for analyses such as collec- 
tive flow and fluctuation which provide important physics at RHIC. Not only 
can pion yields be obtained by statistically fitting of the histogram, pion iden- 
tification can also be obtained track by track using a dE/dx cut on the na^. 
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Fig. 10 illustrates the contamination to pions in percentage from kaons and 
(anti-)protons when we select na n greater than a given value (threshold). We 
can achieve > 95% pion purity at 50% pion efficiency when requiring na^ > 0. 
As have been shown previously, ncrj is actually not centered at zero due to 
insufficient calibration. However, this can be easily corrected for. If one can 
tolerate a pion contamination of 10%, then the efficiency can be as high as 
~ 80%. These represent the worst-case scenario since at higher pt the baryon 
enhancement is less and the contamination should be smaller. 



3. 5 Contamination in proton identification 




1 1.5 2 2.5 3 3.5 4 4.5 5 



p T (GeV/c) 

Fig. 11. m 2 resolution as a function of pr- 

The contamination in (anti-)proton identification is a little more complicated. 
As shown in Fig. 2, the separation of dE/dx between protons and kaons is 
~ la at intermediate/high px- Therefore the purity of protons at p? < 4 
GeV/c is mainly driven by the m 2 resolution of TOFr. The m 2 resolution is 
formulated by: 



2 _ 6p 2 p 2 25t p 2 25L 
dm ~Wl~ 2 ®J 2 ~®J 2 ~' 

where p is the momentum, t is the time of flight and L is the track length of 
the particle. Fig. 11 depicts the m 2 resolution for all charged particles as a 
function of Pt, at —0.5 < rj < 0. The momentum resolution is taken from [3]. 
St and 5L are chosen to be llOps and 0.5cm, respectively. At intermediate/high 
Pt, the m 2 resolution increases approximately quadratic with momentum. If 
we require m 2 > 0.88 (mass square of p(p)) and na p ^ < 0, which means an 
(anti-)proton efficiency of 25%, the contamination from pions and kaons is 
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dE/dx distribution in minbias 0.5<lnl<0.9, 4.0<p_T<4.5 GeV/c 



dE/dx distribution in minbias d+Au 4.0<pT<4.5 GeV/c 
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Fig. 12. dE/dx distribution normal- 
ized by pion dE/dx and offset by 
±5 for positive and negative charge 
at 4 < p T < 4.5 GeV/c and 
0.5 < \t]\ < 0.9, respectively. The 
distribution is from minimum-bias 



Fig. 13. dE/dx distribution normal- 
ized by pion dE/dx (4 < px < 4.5 
GeV/c) and offset by ±5 for inclu- 
sive hadrons at 0.5 < i] < 1.0 and 
— 1.0 < 7] < —0.5, respectively. The 
distribution is from minimum-bias 



Au+Au collisions at y/s NN = 62.4 



d+Au collisions at y/s NN = 200 GeV. 



GeV. 



~ 10% at 4 < pt < 5 GeV/c. We have used the invariant yields of 7r/K/p 
from [7] [13] in this estimation. 

At higher (> 5 GeV/c), TOFr is not effective in identifying (anti-)protons 
anymore. The contamination to (anti-)proton identification from kaons is 
mainly determined by the >lcr separation between them. Therefore, at 50% 
efficiency (na p (p) < 0), kaon contamination is about 20% assuming kaon and 
proton yields are equal. Again we can use different methods to evaluate the 
contamination and efficiency as described in previous section. 

At STAR, the resolution of dE/dx measurement depends on many factors, 
including magnetic field setting, event multiplicity, beam luminosity, track 
length and drift distance. It also depends on the number of hits in TPC 
used to calculate the ionization energy loss for a given track. Excellent read- 
out electronics system (pulse shape control, fine linearity and large dynamic 
range, etc.) and careful offline calibration are important to achieve high-quality 
dE/dx measurement. Particle identification may have difference efficiency and 
purity, depending on the experimental setting and cuts used in analysis. Ba- 
sically, dE/dx resolution is improved with longer track length, shorter drift 
distance, stronger magnetic field, lower multiplicity and beam luminosity, as 
well as more hits in TPC used for track fitting. Fig. 12 and 13 shows pion 
identification at 0.5 < \r}\ < 0.9 in Au+Au collisions at y/s NN = 200 GeV 
and 1.0 > \t]\ > 0.5 in d+Au collisions at 4 < p T < 4.5 GeV/c, respectively. 
Due to longer track and shorter drift distance for particles produced at higher 
|?7|, the dE/dx resolution gets better. Thus the separation between pions and 
Kaons or(anti-)protons in Fig. 12 and Fig. 13 are larger than that in Fig. 3, 
and better identification can be achieved. From the relative heights of pion 
and proton peaks, it is also obvious that the particle compositions in Au+Au 
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at y / s^"=62.4 GeV and d+Au at v /s^"=200 GeV are different. 
4 Electron Identification at low/intermediate p T 




Fig. 14. dE/dx in the TPC vs. particle momentum (p) without (upper panel) and 
with (lower panel) TOFr velocity cut of \l/(3 - 1| < 0.03. 

In addition to its hadron PID capability, TOFr detector can be used to iden- 
tify electrons by combining it with dE/dx from TPC. We were not able to 
measure dilepton with the prototype MRPC TOF tray installed in STAR due 
to small TOF acceptance. However, single electron spectra are sensitive to 
charm production and future azimuthal 2ir coverage of TOF will enable us 
to do dilepton physics[8]. In this section, we will discuss the electron identi- 
fication method, the hadron rejection power at low pr and rejection of the 
photonic background from 7 conversions and Dalitz decays [14]. 

The top panel of Fig. 14 shows the 2-D scatter plot of dE/dx as a function of 
the particle momentum (p) for charged particles with TOFr matched hits from 
d+Au collisions. The bottom panel shows that slow hadrons are eliminated 
and the electron band is well separated from the hadron band with a particle 
velocity (/3) requirement at \l/(3 — 1\ < 0.03. Electrons can then be selected 
with the following cut (1): 

dE/dx{p) > 2.4 + 0.65 x (1 - e -(P-°- 15 )/ ai ) + 0.1 x p (1) 
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dE/dx (keV/cm) dE/dx (keV/cm) 

Fig. 15. dE/dx distribution with TOFr velocity cut in 2 px bins. The arrows de- 
note the cut from Eq.( 1). The curve in left plot (0.4 < px < 0.5GeV/c) shows 
a Gaussian plus exponential fit to the dE/dx distribution, while that in right plot 
(2.0 <pr < 2.5GeV/c) is a 2-Gaussian fit. 

where p is in GeV/c and dE/dx is in keV/cm. With the combination of dE/dx 
from TPC and /3 from TOFr, electrons can be identified above p ~ 0.15 GeV/c, 
while the high pt reach is limited by the statistics in this analysis. 

4-1 contamination for electron identification 

Hadron contamination to electron identification with TOFr velocity cut was 
studied from the dE/dx distribution in each pt bin. A Gaussian function with 
an exponential tail is used in the fit. At px — 2 — 3 GeV/c, 2-Gaussian fit is 
also performed and shows little difference from the Gaussian plus exponential 
fit. Fig. 15 shows the results in two px bins from d+Au collisions. The arrows 
denote the cut from Eq. ( 1). Hadron contamination ratio is estimated from 
these fits, and shown in Fig. 16. The electron efficiency is obtained by vary- 
ing the cut in Eq. ( 1). At low px, high electron efficiency and low hadron 
contamination can be achieved. Even at intermediate px, the contamination 
can be reduced to a low level (~ 1%) if a lower electron efficiency (~ 50%) is 
chosen. Fig. 17 shows electron to hadron ratio vs. px before and after TOFr 
velocity selection. The hadron rejection power can be evaluated by (hadron 
contamination) x(e/h) /(electron efficiency) for a specific condition. It was es- 
timated to be at 10~ 5 level at px < 1-0 GeV/c. At higher p T , additional hadron 
rejection from electromagnetic calorimeter in similar fiducial coverage can help 
reach the same rejection power [16]. 

4-2 Single electrons from Au+Au collisions at y/s NN = 62.4 GeV 

The purpose of the data analysis of d+Au and p+p collisions is to set up 
the baseline for heavy ion collisions. Since the large data sample of Au+Au 
collisions at y/s NN = 200 GeV is not available yet, the smaller sample of Au+Au 
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Fig. 16. The contamination from hadrons to electron identification as a function of 
electron efficiency in different pt range, with TOFr velocity cut. 
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Fig. 17. Electron/hadron ratio vs. transverse momentum (pr), with and without 
TOFr velocity cut. 

collisions at y/s NN = 62.4 GeV has been analyzed to understand the necessary 
techniques for Au+Au collisions. There is a significant increase of multiplicities 
in Au+Au collisions as compared to that in p+p and d+Au collisions, and a 
consequent decrease of track quality. In this section, we discuss the electron 
measurements from TOF and dE/dx method and the photonic background 
estimation. Since the charm yield is low compared to photonic contributions at 
y/s NN = 62.4 GeV, any excess of electrons above the photonic background with 
acceptable errors is not expected. This provides a testing ground to evaluate 
the quality of background reconstruction. 



STAR has accumulated ~ 15 million events in a relatively short run with 
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Au+Au collisions at y/s NN = 62.4 GeV. With the minimum bias trigger (0 — 
80%) and vertex z position selection, the events used in this analysis are ~ 6.4 
million. The resolution for TOF detectors is ~ 110 ps for TOF system, with 
~ 55 ps start timing resolution included. The hadron PID capability was 
reported in [7] . As in d+ Au and p+p collisions, electrons can be identified by 
combining TOF and dE/dx in the TPC. Electrons were selected to originate 
from the primary interaction vertex according to the criteria shown in Table 
2. 

Table 2 

Electron selection criteria in Au+Au 



Method 


TOF+dE/dx 


\VertexZ\ < 


30 cm 


primary track ? 


Yes 


nFitPts > 


25 


ndEdxPts > 


15 


rapidity 


(-1.0, 0) 


X 2 /ndf 


(0., 3.0) 


(3 from TOF 


\l/(3 - 1| < 0.03 


TOFr hit quality 


30 < ADC < 300 




-2.7 < z loca i/cm< 3.4 




\yiocai - yc\ < 1-9 cm 


TOFp hit quality 


thi < ADC < th 2 




2.0 < z local /cm< 18.0 




0.4 < y local /cm< 3.2 



Fig. 18 shows the dE/dx vs. particle momentum after a (3 cut from TOF. 
The electron dE/dx band can be separated from that of hadrons. The dE/dx 
resolution in Au+Au collisions decreases when compared to that in d+Au 
and p+p due to much higher multiplicities and a multiplicity dependence 
of dE/dx calibration which is not currently implemented. We fit the dE/dx 
distribution around electron peak with both two-gaussian function and expo- 
nent ial+gaussian function to extract the electron raw yields in each p T bin. 

Hadron contamination becomes larger at p T > 1-5 GeV/c if we select elec- 
trons with the same efficiency as in p+p and d+Au collisions. Fig. 19 shows 
the hadron contamination fractions under different electron selections. The 
selection a e > with electron efficiency of about 50% is necessary to have the 
contaminations less than 10% at Pt > 1-5 GeV/c. At lower p T , the identifica- 
tion is not affected as much by the high multiplicity since most of the hadrons 
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Fig. 18. dE/dx vs. particle momentum after a TOF (3 cut (|l//3 — 1| < 0.03) 
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Fig. 19. Hadron contamination fractions for different electron selections in Au+Au 
62.4 GeV collisions. 



are rejected by time-of-flight selection. 



4-3 Estimate of electrons from photonic source in Au+Au 



Photon conversions 7 — > e + e~ and 7r° — > ^ye + e~ Dalitz decays are the domi- 
nant photonic sources of electron background. Since only one TOF tray (120th 
of the proposed TOF) is available, electron identification of both e + e~ daugh- 
ters from conversion or Dalitz decays is not possible at this momentum. To 
measure the background photonic electron spectra, the invariant mass and 
opening angle of the e + e~ pairs were constructed from an electron (positron) 
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Table 3 

Partner candidate selection criteria in Au+Au 



charge 


opposite to tagged track 


primary/global ? 


global 


nFitPts > 


15 


nFitPts/nMax > 


0.52 


X 2 /ndf 


(0., 3.0) 


cr e 


(-1., 3.0)t 


dca of e + ,e~ 


(0.0, 3.0) cm 



t several different a e cuts have been studied. 




0.5 1 1.5 2 

Invariant Mass M e+e (GeV/c 2 ) 



Fig. 20. Main photonic background reconstruction in Au+Au collisions. The crosses 
depict the real electron pair candidate invariant mass distributions and the his- 
tograms represent the combinatorial backgrounds. The photon conversion peak is 
clear in ~ mass region, while the ir° Dalitz is hard to see. 

in TOFr and every other positron (electron) candidate reconstructed in the 
TPC [15]. A secondary vertex at the conversion point was not required. In 
Au+Au collisions, due to large multiplicity, the partner track reconstruction 
will lead to a large combinatorial background due to hadron contamination 
without an additional PID from TOF. Even with a stringent dE/dx selection, 
this combinatorial background is still significant. In the future, dilepton stud- 
ies would be possible with 2tt coverage of TOF for selecting clean e + e~ pairs. 
Fig. 20 shows the electron pair candidate invariant mass distribution. The 
tagged electron was selected from TOF with the cuts shown in Table 2 and 
additional < na e < 3. The partner track candidate was selected according 
to Table 3. 
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Fig. 22. The reconstruction effi- 
ciency-corrected photonic back- 
ground over the raw inclusive 
electron yield under different elec- 
tron/positron track selections. The 
numbers in the brackets on the plot 
show the a e cut for the tagged track 
and the partner track, respectively. 



Fig. 21. Photonic background recon- 
struction efficiency from Au+Au 62.4 
GeV HIJING simulations. Also shown 
on the plot is that from d+Au 200 
GeV HIJING simulations. 



The combinatorial background distribution was provided by rotating the part- 
ner track momentum J? — > — jf, and normalized to the distribution of the real 
electron pair candidates in the region 0.8 < M e + e -/< 2.0 (GeV/c 2 ) in the in- 
variant mass spectrum. Fig. 20 shows the results in 0.4 < px < 0.6 GeV/c 
where px is the px of tagged electrons. The plot shows that the combinatorial 
background was very well reproduced. The peak from photon conversion in 
detector material is clearly visible near zero mass region, and its offset from 
zero is due to the opening angle resolution [15] in TPC tracking. The n° Dalitz 
contribution is not visible in this case, possibly because the Dalitz contribu- 
tion is much smaller and distribution is much broader compared to conversion 
processes. We subtracted the combinatorial background from the real distribu- 
tion, and integrated the remaining distribution from — 0.15 GeV/c 2 to get the 
reconstructed main photonic background raw yield. In this case, we assumed 
that both photon conversion and 7r° Dalitz decays were reconstructed. 

The background raw yield need to be corrected for the reconstruction ef- 
ficiency, as we did it for d+Au and p+p collisions [8]. This efficiency was 
calculated from Au+Au 62.4 GeV HIJING events with full detector MC sim- 
ulations. After \Vz\ < 30 cm cut, ~ 53 K events were used in the calculation. 
We took all TPC electron tracks without any dE/dx and TOF hit cut to im- 
prove the statistics. The procedure is the same as we did in p+p and d+Au 
collisions [8]. Fig. 21 shows the background reconstruction efficiency in Au+Au 
62.4 GeV compared with d+Au results. Due to the higher multiplicities, the 
reconstruction efficiency is relatively lower in Au+Au than that in d+Au. 

This efficiency was used to correct for the photonic background raw yield ob- 
tained from above. In addition, a ~ 5% fraction of other photonic background 
(from d+Au results) and the na e selection efficiency were also included. The 
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corrected background spectrum is compared with the inclusive spectrum as 
shown Fig. 22. The reconstructed background matches the total inclusive elec- 
tron spectrum. This is expected since the charm yield is more than an order of 
magnitude lower than the inclusive electron yield at 62.4 GeV [14]. This also 
demonstrates that we can reconstruct fully the photonic background from 7 
conversions and Dalitz decays. 

The performance of the electron identification and background studies for 
Au+Au system makes us confident on extracting the charm signal at 2 — 3 
GeV/c from the coming 30 M minimum bias 200 GeV Au+Au data. With the 
STAR upgrade of inner tracker -Heavy Flavor Tracker (HFT), we will be able 
to reject photon conversions by requiring additional hits in HFT, to reject 
electrons from charm semileptonic decays by its displaced secondary vertex, 
and to further reconstruct /reject the low Pt partner of the electrons from 
the Dalitz decays using inner tracker only. Electron identification at px > 1 
GeV/c can be improved by combining dE/dx from TPC, velocity from TOF 
and energy from electromagnetic calorimeter. Detailed studies are underway 
and are beyond the scope of this article. 



5 Conclusion 

In summary, we have developed a technique to extend the particle identi- 
fication up to high pt at STAR. By combining information from the TPC 
and TOF, we can measure pion, kaon and (anti-)proton in the intermedi- 
ate/high p T range and electron in low/intermediate px (high p T limited by 
statistics). Preliminary spectra obtained in Au+Au collisions at v /s^~=62.4 
GeV shows reliable particle identification to pt ~ 7 — 8 GeV/c for pions and 
(anti-)protons, and ~ 3 GeV/c for kaons, respectively. The electrons were 
identified at 0.15 < pr < 4 GeV/c, in d+Au, p+p collisions at v /s^J7=200 
GeV and Au+Au collisions at v /s^7=62.4 GeV. A purity of over 95% for pion 
identification and ~ 70 — 90% for proton identification was obtained with this 
method at intermediate/high p?- We achieved a hadron rejection power at 
10~ 5 level for the electron identification at low pt- 
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